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This article describes a digital Doppler extraction demonstration with the Ad- 
vanced Receiver II (ARX II) tracking Pioneer 10 and Voyager 2 . The measured 
results are compared with those of the Block IV receiver that was operating in 
parallel with the ARX IL It is shown that the ARX II outperforms the Block IV 
receiver in terms of Allan variance of the Doppler residuals, the amount of which 
depends on the scenario of interest. 


I. Introduction 

The Advanced Receiver (ARX) [1] was modified and 
upgraded to accommodate higher telemetry rates and to 
enable digital Doppler extraction. The previous receiver 
(ARX I) underwent several telemetry demonstrations [2-4] 
at Goldstone, but was unable to perform Doppler extrac- 
tion because the sampling clock was not fixed, but was an 
integer multiple of the symbol rate. This is referred to as 
“synchronous sampling” (synchronous with respect to the 
symbol epoch). In the new receiver (ARX II) [5], the sam- 
pling clock is fixed and is derived from the Frequency and 
Timing Subsystem (FTS). Hence, sampling is performed 
asynchronously with the symbol period and may result in 
a noninteger number of samples per symbol. However, a 
fixed time base is present in the system and digital Doppler 
extraction can be performed. 


The architecture of the ARX II is discussed in [5] 
along with the various signal-processing schemes employed 
to perform the functions of carrier and subcarrier tracking 
and symbol synchronization, among others. In this article, 
only the phase-locked loop is described, as it was used to 
track the residual component of the received signal from 
the two spacecraft. 

Both analog and digital loop closures can be mecha- 
nized with the ARX II and both are depicted in Fig. 1. 
A software command can switch between both methods of 
loop closure without losing the carrier lock. 

In the analog loop closure, the incoming signal is at 
70 MHz with about 17-MHz intermediate frequency (IF) 
bandwidth. It is mixed with an 80-MHz locked reference 
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to produce the 10-MHz IF, which is sampled at 40 MHz us- 
ing an 8-bit an alog-to- digital converter (A/D). The down- 
conversion to baseband in-phase and quadrature (I and 
Q) samples is then accomplished by digital mixing using 
look-up tables. The baseband signal has about an 8-MHz 
baseband bandwidth and thus a decimation by two is per- 
formed to process the signal at the sufficient 20-MHz pro- 
cessing rate (effective sampling rate). The Q samples are 
accumulated to reduce the rate to the more appropriate 
loop update rate (50 Hz-1 kHz) to enable a software imple- 
mentation of the loop filter F(z). The filter output is used 
to adjust the frequency of a numerically controlled oscilla- 
tor (NCO), which nominally operates at 4 MHz. The out- 
put of the NCO drives a digital-to-analog converter (DAC), 
producing a 4-MHz analog sine wave, which is mixed with 
an 84-MHz fixed reference to produce the 80-MHz locked 
signal that closes the loop. 

The operation of the digital loop closure is similar 
to that of the analog loop closure, except that the 80-MHz 
signal is a fixed reference and an NCO (running at 10 MHz) 
is used to close the loop. In this case, the A/D operates on 
a signal that is not centered at exactly 10 MHz, but rather 
shifted according to the instantaneous Doppler. 

The ARX II accepts an input signal in the 200-400- 
MHz range from the IF Distribution Assembly output and 
downconverts it to the 70-MHz range using an external 
synthesizer (typically a Marconi), as shown in Fig. 2 for the 
digital loop closure. A total-power automatic gain control 
(AGC) circuit operates on the signal to ensure that the 
resulting 10-MHz IF signal lies within the dynamic range 
of the A/D. 

The residual component of the received signal at the 
output of the AGC can be modeled as 

r{t) = y/2F c sm(uit + e c ) + n(t) ( 1 ) 


at 40 MHz. Following I and Q mixing, the samples are 
accumulated to produce 

/ cr (n) = \/^cos<£ c (n) + nj(n) (3) 

Qcr(n) = y/K sin <t>c(n) + «o( n ) ( 4 ) 


where n/(n) and HQ(n) are independent Gaussian ran- 
dom variables with variances No/2T Ui where T u denotes 
the loop update period and <f> c (n) the carrier phase error. 
The Q cr samples update the loop filter, which is given by 


F(z) = G\ + 


G 2 ° 3 

1 - z- 1 (l-z- 1 )* 


(5) 


where 


e£ 

11 

(6) 

G 2 = rd 2 /T u 

(7) 

G 3 = krd a /T u 

(8) 

d = AB l T u (r - i)/r(r - k + 1) 

(9) 


where Bl is the design loop bandwidth in hertz, r is typ- 
ically 2 or 4 and is equal to 4£ (where £ is the damping 
ratio), and k is a type-III loop gain parameter (k = 0 for a 
type-II loop) with typical values ranging from 1/4 to 1/2. 

Digital Doppler extraction consists of recording the 
phase estimate of the incoming waveform time tagged to 
station time using a Time Code Translator (TCT). The 
resolution of Doppler data depends heavily on the number 
of bits used in the NCO to represent phase and frequency. 


where P c denotes the carrier power, the IF frequency 
(70-MHz) in rad/sec, and B c the carrier phase in radians. 
The data modulation is ignored here as only the residual 
carrier component is of interest. The narrowband noise 
n(t) can be written as 

n(t) = y/2n c (t) co&(u>it -f 0 C ) - V2n s {i) sin(w,< + 0 C ) (2) 

where n c (i) and n a (t) are statistically independent, sta- 
tionary, band-limited, white-Gaussian-noise processes with 
one-sided spectral density N 0 watts/Hz and one-sided 
bandwidth W Hz (roughly 17 MHz in this case). The 70- 
MHz IF signal is downconverted to 10 MHz and sampled 


In the digital loop closure, the NCO accumulates frac- 
tional phase using 32 bits, but only the most significant 12 
bits are used in the sin/cos look-up table. This results 
in a phase resolution of 0.24 mcycle (2“ 12 ). The phase is 
recorded, external to the NCO, using 48 bits for the in- 
teger part and 12 bits for the fractional part. The NCO, 
in this case, operates at a nominal frequency of 10 MHz 
and is clocked by a 40-MHz clock. This results in 9.3-mIIz 
(4 x 10 7 /2 32 ) frequency resolution. 

In the analog loop closure case, the fractional phase is 
recorded using 14 bits, two of which are accumulated exter- 
nal to the NCO. The resulting phase resolution is 0.06 mcy- 
cle, better than that of the digital loop closure. The NCO, 
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in this case, is actually operating at a nominal frequency of 
16 MHz with a 20-MHz clock, but it effectively produces a 
4-MHz signal because of the additional bits. The resulting 
frequency resolution is 1.16 mHz (2 x 10 7 /2 3 ^), better than 
the digital loop case. 

However, the improved resolution for the analog loop 
case is at the expense of limiting the maximum nominal fre- 
quency to one-fourth of the clock rate (the practical limit 
of an NCO is roughly 40 percent of the clock rate). That 
in turn limits the dynamic excursion of the frequency to 
±1 MHz in the analog loop closure, whereas it is about 
±5 MHz in the digital closure. Time tagging is performed 
at ±1 /isec accuracy with a ±1 nsec stability. The logging 
rate can be chosen by software from 1, 2, 5, 10, or 20 Hz. 

II. Demonstration Description and Results 

The demonstration consists of operating the ARX II 
in parallel with the Block IV receiver; both track the same 
spacecraft. Both receivers record data for an identical pe- 
riod of time. The data are later analyzed to quantify and 
compare the respective performances (see Fig. 3). In all 
the results to follow, the ARX II was operating with a 
500-Hz loop update rate. 

A. Allan Variance 

The X-band downlink of Voyager 2 (channel 14, 
8415.000000 MHz) and the S-band downlink of Pioneer 10 
(channel 7, 2292.407407 MHz) were tracked to determine 
the relative effectiveness of the ARX II and the Block IV 
receiver in Doppler measurements. 

Doppler residuals from four data arcs, as defined in 
Table 1, are reported here. The first arc is from the 
Voyager-2 X-band downlink tracked on Greenwich Day of 
Year (DOY) 195, the second from the Pioneer-10 downlink 
tracked on DOY 195, and the last two from the Pioneer-10 
downlink tracked on DOY 201. During all of these tracks, 
the spacecraft operated in a three-way Doppler mode, thus 
eliminating problems associated with the spacecrafts’ on- 
board oscillators. For each of these arcs, the downlink was 
simultaneously tracked by the ARX II and a Block IV re- 
ceiver. The downconverted carrier phase as measured by 
the ARX II was recorded once per second and the Doppler 
count as measured by the Block IV Doppler extractor was 
recorded once per second. 

Use was made of two time scales during the analysis: 
a 24-hour clock for Greenwich Mean Time and the number 
of seconds past the January 1, 1950 epoch. This “J 50” 
time was needed to do a linear interpolation on the uplink 


frequencies as reported in the Program Frequency Records 
(Table 2). 

The Doppler observable is calculated for the ARX IT 
as 

Ido = 48 G/e + Jnco — /i — fi — h (10) 

where foo is the Doppler observable, G the transponder 
ratio (880/221 for the Voyager track and 240/221 for the 
Pioneer tracks), / e the exciter frequency, fsco the differ- 
ence in phase as measured one second apart by the ARX II, 
fi the frequency of the first local oscillator (8.1 GHz for the 
Voyager track, 2.0 GHz for the Pioneer tracks), fa the fre- 
quency of the second local oscillator, and finally, / 3 , which 
is 80 MHz for digital closure and 74 MHz for analog closure 
of the ARX II loop. 

The exciter operates at approximately 44 MHz, 
whereas the Program Frequency Records report the exciter 
frequency divided by two (approximately 22 MHz). The 
exciter frequency is found by interpolating the frequency 
ramps of Table 2 and then multiplying by two. The sec- 
ond local oscillator is a Marconi synthesizer; its setting for 
each of the tracks is given in Table 3. The final term in 
Eq. (10) depends on whether the ARX II had digital or 
analog loop closure; this information is given in Table 4. 

Prior to DOY 201, there was a hardware clocking 
problem that caused the cycle count recorded by the 
ARX II to be low by exactly one cycle on some (randomly 
occurring) samples and these were processed by a cycle 
correction algorithm that was used to correct the data. 
The DOY 201 data did not need to be processed by this 
algorithm because the hardware problem was fixed by that 
date. 

The Doppler residual is calculated from the Doppler 
observable as 

/a = Ido - hi ~ It (11) 

where /a is the Doppler residual, / m the modeled (based 
on trajectory) Doppler observable, and fx the tropospheric 
correction. 

Equation (11) applies to both the ARX II Doppler 
residual and to the Block IV Doppler residual. The 
Block IV Doppler observable, the modeled Doppler ob- 
servable, and the tropospheric correction are all available 
on the REG RES file provided by navigation personnel. 
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Doppler residuals, spaced one second apart, were cal- 
culated for both receivers for each of the four data arcs 
defined in Table 1. The Doppler residual record for part 
of the Voyager 195 data arc is shown in Fig. 4. In that fig- 
ure, a bias was added to the Block IV residuals to displace 
them from the ARX II residuals so that the two traces 
could be discerned. In fact, the means of those two traces 
are in excellent agreement. The mean values of the Doppler 
residuals for the ARX II and for the Block IV receiver are 
given in Table 5. It is noted that the two receivers agree 
on the magnitude of the Doppler effect to less than 1 mHz. 

The square-root Allan variance is used here to char- 
acterize the jitter in the Doppler residual measurements. 
It is defined as 

( Ty(r ) = -7=- \J < [fb(t + T ) / a (*)] 2 > ( 12 ) 

v 2v 

where <Ty(r) is the square-root Allan variance, t the inte- 
gration time, v the downlink frequency, and < • > denotes 
the time average over the entire data arc. The square- 
root Allan variances computed for a one-second integra- 
tion time are contrasted in Table 6. Also given in Table 6 
are modeled values of that quantity based on the assump- 
tion that downlink thermal noise dominates the frequency 
instability for such a short integration time. The modeled 
values came from [6], i.e., 


<ry( r ) = 


1 

2wtv 



(13) 


where P c /N 0 denotes the carrier power-to-noise spectral 
density ratio (from Table 4) and Bi the one-sided loop 
bandwidth. 

The modeled and measured variances for short inte- 
gration times as presented in Table 6 agree closely for the 
Pioneer tracks, indicating that the Doppler measurement 
error really is dominated by downlink thermal noise. For 
the Voyager 195 track, the measured variance is somewhat 
higher than the modeled variance; this could be due to the 
carrier power not being constant during the pass. 

The square-root Allan variance of the Doppler resid- 
ual is plotted as a function of integration time in Figs. 5 
and 6 for Voyager and Pioneer, respectively. These plots 
show a thermal noise dominance with reciprocal r depen- 
dence for the smaller integration times, and the advantage 
of the ARX IPs smaller loop bandwidth is made plain. For 
larger integration times, other error sources dominate and 
the bandwidth advantage of the ARX II is lost. With the 


Pioneer data, the two traces actually begin to coalesce, 
indicating a common error source for the two receivers. 
This common error source, which dominates for larger in- 
tegration times, is almost certainly solar plasma. The Sun- 
Earth-probe angle was 38.2 deg for Pioneer on DOY 195. 
The Pioneer 201 data arcs are not long enough to pro- 
vide a statistically significant Allan variance for the longer 
integration times. 

B. Receiver Tracking Performance 

The ARX II acquires the carrier using an automated 
procedure, which consists of performing several software 
fast Fourier transforms (FFTs) around the frequency pre- 
dict. The size and the rate of the FFTs are controlled by 
software, and the spectra are noncoherent^ averaged to 
reduce the noise. The largest peak is then selected and the 
loop enabled. The acquisition time depends on the speed 
of the FFT computation. A future upgrade is planned that 
will use a dedicated FFT board to perform the computa- 
tions. This will dramatically reduce the acquisition time 
to less than a second. 

For performance evaluation, the ARX II records sev- 
eral parameters. These include: ^carrier power-to-noise, 
spectral, density-level estimates, P c /N 0 ] averaged and nor- 
malized I and Q samples; steady-state phase-error esti- 
mates <f> SS) etc. 

The averaged I sample can be used in real time to 
monitor the loop gain and, hence, the loop bandwidth. 
With proper normalization, the averaged I samples should 
be unity, assuming that the amplitude of the incoming 
signal is “close” to its predict and, hence, that the operat- 
ing loop bandwidth is its designed value. That averaging 
should be performed when the steady-state error due to dy- 
namics is negligible; otherwise, no information about loop 
gain can be deduced. A steady-state phase-error estimate 
can be obtained in real time by computing 


= Ian 


/ E^ :+£+1 ar(n) \ 

VeS +£+ 1 M")/ 


(14) 


where L controls the estimation period. In the ARX II, 
£(j) can be sent to the monitor’s screen at different rates 
ranging from once every second to once every minute. An- 
other key parameter is the incoming Pc/Nq } which is also 
estimated by the receiver according to 



(En=r L+ 1 Mn )) 2 
2r u En=j + L + 1( 5?r(n) 


(15) 


163 


Figure 7 depicts the various parameters for a Pioneer 
track using analog loop closure with a type-III loop. In this 
case, the estimation period was 10 sec and the designed 
loop bandwidth about 0.6 Hz. It is clear that the aver- 
aged I is close to unity and the averaged Q close to zero, 
which indicates that the operating bandwidth is 0.6 Hz, 
as designed. The estimated P c /No is roughly 12.5 dB-Hz 
(compared with a 13 dB-Hz predict) and the steady-state 
phase error is approximately zero. 

In Fig. 8, Pioneer results are also shown for the digital 
loop closure case with a type-II loop. The averaged Q is 
no longer zero as expected, due to the 7-deg steady-state 
phase error (since type-II loops cannot track Doppler rates 
with zero phase error). The resulting P c /Nq estimates are 
also degraded by about 0.6 dB to yield roughly 11.9 dB- 
Hz. Finally, Fig. 9 depicts the results for a Voyager 2 
track using a type-III loop with digital closure. Note that 


the normalized I is not unity, which implies that the op- 
erating loop bandwidth is actually 2.5 Hz (1.2 times the 
designed bandwidth of 2 Hz). The steady-state error is 
zero as expected and the estimated P c /N o is about 31.3 
dB-Hz (compared with a 32-dB-Hz predict). The estima- 
tion period for the Voyager data is one second due to the 
relatively “strong” signal. 

III. Conclusion 

This article described digital Doppler extraction with 
the ARX II and the results were compared to the Block IV 
receiver Doppler extractor. It was shown that digital 
Doppler extraction did not result in any degradation rela- 
tive to Block IV, and in most cases performed better with 
the ARX II. From the standpoint of Doppler extraction 
performance, ARX II is an excellent candidate for future 
Block V receiver implementation. 
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Table 1. Definition of tracks 


Track 

Start time 

End time 

hh:mm:ss 

J 50 Time, sec 

hh:mm:ss 

J 50 Time, sec 

Voyager 195 

06:41:00 

1247553660 

07:52:00 

1247557920 

Pioneer 195 

14:11:00 

1247580660 

14:52:00 

1247583120 

Pioneer 201 A 

12:18:00 

1248092280 

12:35:00 

1248093300 

Pioneer 201 B 

14:46:00 

1248101160 

15:07:00 

1248102420 


Table 2. Data from Program Frequency Records 


Pass Time, J 50, sec Start frequency, Hz Ramp rate, Hz/sec 


V195 1247552633 
V195 1247557733 
P195 1247573336 
P201 1248089640 


0,22015336732360 SOD + 08 
0.2201 534843686050D + 08 
0.21 984 78999999950D + 08 
0.21 984 6 29999999 SOD + 08 


0.2295000000000000D - 02 
0.2235000000000000D - 02 
0. 0000000000000000 D + 00 
0.00000000000 OOOOOD + 00 


Table 3. Second local oscillator frequency 


Track 


Second local oscillator (.fo), Hz 


Voyager 195 
Pioneer 195 
Pioneer 201 A 
Pioneer 20 IB 


244416020 

222078470 

222094640 

222094640 



Table 4. Receiver parameters 


Track 

P c /No, dB-Hz 

Block IV 
One-sided loop 
bandwidth, Hz 

ARX II 
One-sided loop 
bandwidth, Hz 

ARX II 
Loop 
closure 

Voyager 195 

30 

54 

2.0 

Digital 

Pioneer 195 

11 

2.5 

2.0 

Digital 

Pioneer 201 A 

11 

2.5 

1.0 

Digital 

Pioneer 201 B 

11 

2.5 

0.6 

Analog 


Table 5. Mean of Doppler residual 



Mean of Doppler residual, Hz 

Track 

Block IV 

ARX II 

Voyager 195 

-0.130 

-0.130 

Pioneer 195 

-0.153 

-0.153 

Pioneer 201 A 

-0.148 

-0.147 

Pioneer 201 B 

-0.149 

-0.149 


Table 6. Square-root Allan variance of Doppler residuals 


1-second (square-root) Allan 

variance of Doppler residuals 

Track 

Block IV 

arx n 


Modeled 

Measured 

Modeled 

Measured 

Voyager 195 
Pioneer 195 
Pioneer 201A 
Pioneer 201 B 

7.6 X 10 -12 
5.4 X 10-" 
5.4 X io- u 
5.4 X 10 -n 

1.1 x 10 -n 

5.1 x 10 -11 
5.4 X 10 -11 
5.4 X 10“" 

1.5 X 10~ 12 
4.8 X lO -11 
3.4 x 10- 11 

2.6 X 10~ n 

2.3 X 10“ 12 
4.7 X 10“ n 

2.9 X lO" 11 

1.9 X 10“ 13 




Fig. 1. Analog versus digital loop closure. 
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S: 2. 2-2. 3 GHz 

X: 8.4-S.5 GHz (8.2-8.6 GHz FOR VLBI) 



Fig. 2. ARX 11 phase-locked loop. 


















FROM MASER 



Fig. 3. Experiment description. 



Fig. 4. Doppler residual for Voyager 2, DOY 195. 














SQUARE-ROOT ALLAN VARIANCE OF DOPPLER RESIDUAL 



INTEGRATION TIME, sec 


Fig. 5. Allan variance of Doppler residual for Voyager 2, DOY 195. 



INTEGRATION TIME, sec 


Fig. 6. Allan variance of Doppler residual for Pioneer 10, DOY 195. 
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